00 



Determination of the Mott insulating transition by 
the multi-reference density functional theory 



o 
o 

^ , K. Kusakabe 



Graduate School of Engineering Science, Osaka University, 1-3 Machikaneyama-cho, 
Toyonaka, Osaka, 560-8531, Japan 



Abstract. It is shown that a momentum-boost technique apphed to the extended 
Kohn-Sham scheme enables the computational determination of the Mott insulating 
^ ' transition. Self-consistent solutions are given for correlated electron systems by the 

first-principles calculation defined by the multi-reference density functional theory, 
' in which the effective short-range interaction can be determined by the fluctuation 

reference method. An extension of the Harriman construction is made for the twisted 
boundary condition in order to define the momentum-boost technique in the first- 
principles manner. For an effectively half-filled-band system, the momentum-boost 
, method tells that the period of a metallic ground state by the LDA calculation is 

' shortened to the least period of the insulating phase, indicating occurrence of the 

^ ■ Mott insulating transition. 

> 

"Tj- ' 1. Introduction 

o ■ 

^ ■ Detection of the Mott insulating transition is a desirable function for the first- 

■ principles calculations, which has been demanded for years. A recently developed 
^ . multi-reference density functional theory with the fluctuation reference method defines 
O ■ a self-consistent first-principles calculation, in which a short-range correlation effect 

is explicitly included. [H [21 [3] This technique is a generalization of the Kohn-Sham 
^ . scheme of the electronic structure calculation. [H [5] Incorporation of the effective many- 

■ body system to determine the total energy and the single-particle charge density of the 
electronic state became possible through 1) introduction of a fluctuation-counting term 
and 2) the reformulation of the exchange correlation energy functional as a residual 
exchange correlation energy functional. The effective Hamiltonian appearing in the 
theory is a kind of the Hubbard model[6] or the Anderson model[7]. The LDA-I-U 
Hamiltonian [8] is derived as an approximation for the residual exchange correlation 
functional. f9] 

In this short paper, I show a technique to determine the Mott insulating phase in 
this first-principles calculation. The method is an application of the momentum boost 
technique known in the literature. [IDl [H] To introduce the momentum boost method in 
the density functional theory, one needs to show the N represent ability. This is done by 
extending the Harriman construction [121 [13] fo^' the twisted boundary condition. A test 
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calculation will be shown using a simple artificial system, which may be represented by 
a one- dimensional Hubbard model. In the last part, I will summarize the present work. 

2. Momentum boost technique 

I consider a Born-von-Karman boundary condition with a twist. Consider an array of 
atoms in one dimension and let L the number of atoms in the direction, which is called 
the X direction below. Formally using a phase 27r$/$o with < $ < $o aiid a unit flux 
$0) I introduce a twisted boundary condition in this x direction. 

^(xi, ■ ■ ■ , (r^ + Lae^, (Xi), ■ ■ ■ , xn) 

= exp (27ri$/$o) ^(a^i, ■ ■ ■ , (r^, a^), ■ ■ ■ , x^v), (1) 
c!v+i,. = exp(27rz<l>/<l>o)cl,,. (2) 

Here, I used a combined coordinate Xj = {rj,aj) with the space and spin coordinates. 

e^; is a unit vector and a is the lattice constant in the x direction. \l/ is a many-body 

wavefunction and is a creation operator defined with a properly determined Wannier 

basis 0i(r). Indeed, determining a Wannier transformation fixing a gauge in the unitary 

transformation, [ll] the set of 0j(r) is uniquely determined in each self-consistent step 

of the extended Kohn-Sham scheme. [2] 

Using another gauge transformation, the twisted boundary condition is represented 

27im 27r$ 

by shifted k vectors, k = ( 1 , ky, kz), for the single-particle Bloch orbitals. 

La La$o 

Here m is an integer within a range of < m < L. From the single-particle energy 
£:„(k), one may construct a tight-binding model written in c]^ and Cj^o-- 

Determining the eigen energy of the effective Fermion system, one has an energy- 
flow diagram as a function of If the system is a metal, the flow given by adiabatic 
connection of the ground state should show a long extended AB period, while the period 
has to be $0; when there is a charge gap due to formation of a Mott gap. Although the 
density-functional theory utilized in the present work is the ground-state formulation, 
one can have a signal of the change in the period as disappearance of a cusp in the 
lowest-energy flow. In case of the Mott insulator, the flux line becomes a smooth curve 
without any cusp. The appearance of the isolated lowest flow lie implies uniqueness 
of the ground state separated by a gap from charge excitations in the whole range 
of twist and thus indicates the formation of the Mott gap. If one has a difficulty in 
finding the change of the period with a three-dimensional sample, one may consider a 
one-dimensional setup using a needle like sample. 

Before discussing the result, I show that the representability of the present 
problem is guaranteed. In the present setup, only the boundary condition in the x 
direction is a twisted one. Thus I consider a slice of the charge density p(x) = p(r) 
fixing y and z coordinate for simplicity. We are seeking for ipk{x) satisfying 

p(x) = 5^A,|^(x)p, (3) 
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where k is an integer, < < 1 and = A^. In the twisted boundary condition, 

k 

however, p(x + La) obeys the periodic boundary condition as 

p{x + La) 

= J daidx2 ■ ■ ■ dxN'^{{ri + Lae^, ai), X2, — - , ^n)* 

X *((ri + Lae^, ai),X2, ■ ■ ■ , xn) 
= J dcridx2 - ■ ■dxN'^{{ri,cri),X2, - ■ ■ ,XN)*'^i{ri,cri),X2, - ■ ■ ,Xn) 
= p{x). (4) 

Here, integration with respect to aj should be interpreted as a summation. Thus, we 
can readily prepare the orbital (pk{x) obeying the twisted boundary condition as, 

Mx) = ^ exp(^A;/(a;) + z27r$/$o), (5) 



/vr 



with 

fix) = ^ I p{x')dx'. (6) 
We can immediately show that a set of (pk{x) forms complete orthonormal and that 

k k 

Using the orbital wavefunction ipk{x) obeying the twisted boundary condition, we can 
show existence of a many-body state |\E') whose coordinate expression is the single Slater 
determinant made on ipk{x). 



3. One-dimensional Hydrogen array 

As a test calculation, I consider a one-dimensional Hydrogen array. The system is 
denoted by an outer unit cell with ten atoms for the many-body calculation and an 
inner unit cell with a single atom for the single-particle problem. The lattice parameters 
for the inner cell are a = 1A and h = c = lOA. The determination of V may be 
achieved by setting a reference calculation. [2j Here, to show the change in the period 
explicitly, I consider [//t as a parameter and do only the self-consistent calculation for 
the extended-Kohn-Sham system. Here, t is the value of the nearest-neighbor transfer 
parameter. Note that the present system is represented by a tight-binding model with 
long-range hopping terms. For the exchange-correlation energy functional, I utilized the 
Perdew-Zunger parameterization of the Ceperley-Alder diffusion Monte-Carlo data. |15] 
The TrouUier-Martins soft pseudopotential is used with the cut-off energy of 20 [Ry] . [T6] 
This setup is confirmed to be enough accurate for the discussion below by increasing 
the parameters. The numerical diagonalization with the Lanczos algorithm is used to 
obtain the many-body state of the first-principles Hubbard model. 
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Figure 1. The ground state energy of a Hydrogen array system as a function of the 
phase $ in the boundary condition. The system has 10 atoms. The value of [//t is 
set to be 0, 1, 2, • • • , 5 with t is the nearest-neighbor transfer parameter. For the case 
of ?7 = (LDA calculation), crossing of the flow lines occur at $ = <&o/2, while the 
crossing becomes anti-crossing due to the charge-gap formation for finite U. 



The result of the momentum boost is depicted in Fig. [31 For [/ = 0, the LDA 
calculation shows a crossing in energy flow lines. In this case, the lowest energy flow 
can be traced, since the constrained LDA calculation fixing the filling of each k point 
is available. Once a finite U is introduced, the many-body calculation automatically 
concludes the lowest branch of the energy flows. In this system, we see a continuous 
change of the energy flow which has an energy gap structure at $ = $o/2. This result is 
qualitatively the same as the single-band Hubbard model with only the nearest neighbor 
hopping. [11] If the fluctuation reference method is precisely applied, a finite value of U 
is expected, since the inter-atomic distance of 2A is in a strong-correlation regime for 
the Hydrogen molecule. [2] Thus the present result gives a concrete test for the method 
of determination of the Mott insulator from the first-principles. 

4. Summary 

We have shown that the momentum boost method is formulated rigorously in the 
density functional theory. The N representability is shown for the twisted boundary 
condition. Using the multi-reference density functional theory (MR-DFT), we are able 
to construct an effective interacting Fermion system, which may undergo the Mott 
insulating transition. The momentum boost technique is applicable for this problem to 
detect the transition. In a realistic system, the formation of the Mott gap by the applied 
pressure or by the effective carrier doping would be seen as a change in the period of 
the lowest energy flow. 
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